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Low-energy doublet neutron-deuteron scattering is described on the basis of the triton bound 

q \ and virtual state parameters — the energies and the nuclear vertex constants of these states. 

(N ■ 

The van Oers-Seagrave formula is derived from the Bargmann representation of the S matrix 
■ for a system having two states. The presence of a pole in this formula is shown to be a direct 

corollary of the existence of a low-energy triton virtual state. Simple explicit expressions 
for the nd scattering length and for the pole of the function k cot 5 are obtained in terms of 
^ ■ the triton bound and virtual state parameters. Numerical calculations of the nd low-energy 

^ ! scattering parameters show their high sensitivity to variations in the asymptotic normalization 

O ! constant of the virtual state C% . The C% value fitted in our model to the experimental result 

O ' for the nd scattering length is C% = 0.0592. 
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1. INTRODUCTION 



The low-energy characteristics of three-body systems were studied previously [1, 2] on the 



^ . basis of the two-body model with the Hulthen potential; in particular, a correlation between 

h : 

the binding energy of three hadrons and the hadron-deuteron scattering length was analyzed 
over a wide region of three-body parameters. An important role of the virtual triton state was 
revealed in studying this dependence (Phillips line) in the region of the experimental values of 
the triton binding energy Et and the doublet nd scattering length 2 a n d- The characteristics 
associated with the T — > d + n decay vertex were also calculated for the bound (T) and virtual 
(v) triton states. The results of the calculations for the position of the virtual triton level, 
B v , as reckoned from the threshold of the elastic nd channel and for nuclear vertex constants 
G\ and G\ characterizing the ground and the virtual triton state, respectively, agree with the 
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relevant results of three-body calculations and with experimental data. The results from [1, 2] 
evince a clear-cut correlation between the low-energy characteristics of two- and three-hadron 
systems. 

It was established in [3-6] that the effective-range expansion of the doublet S-wave phase 
shift for nd scattering involves a pole situated in the nonphysical region near the threshold. 
For the function k cot 5, the corresponding four-parameter representation involving a pole is 
referred to as the van Oers-Seagrave formula. It should be noted that this formula was found 
in a purely empirical way without any theoretical justification. 

In order to analyze two- and three-hadron systems, we propose here an approach based 
on the S matrix corresponding to so-called Bargmann potentials [7, 8] . We will consider the 
relation between the ^-matrix representation introduced by Bargmann [7, 8] and the effective- 
range expansion and demonstrate that, on the basis of this representation, the effective-range 
approximation can be obtained for the case where the system in question has one state, bound 
or virtual. An example of such a situation is provided by neutron-proton scattering. We will 
then generalize this analysis to systems having two states, as in the case of doublet neutron- 
deuteron scattering. It will be shown that the van Oers-Seagrave formula for the function 
k cot S involving a pole follows directly from the Bargmann representation of the S matrix for 
a system having two states. This analysis makes it possible to relate the parameters of the 
effective-range expansion to the characteristics of the bound and virtual states of the system. 

We assume that any state of the system is specified by two parameters, the energy corre- 
sponding to the pole of the S matrix on the imaginary axis in the complex plane of the wave 
number k and the nuclear vertex constant, which is directly expressed in terms of the residue 
of the S matrix at this pole [9] . Nuclear vertex constants are fundamental physical character- 
istics of nuclei like more conventional quantities, including mass, spin, and parity. Kinematical 
factors apart, nuclear vertex constants are related to the on-shell amplitude for the virtual 
or real decay (or fusion) of a nucleus into two fragments. The properties of vertex constants, 
their values for a number of nuclei, and methods for determining them experimentally and 
theoretically were surveyed elsewhere [9]. 
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2. DERIVATION OF THE EFFECTIVE-RANGE APPROXIMATION FROM 
THE BARGMANN REPRESENTATION OF THE S MATRIX INVOLVING 
ONE STATE: NEUTRON-PROTON SYSTEM 



Bargmann [7] proposed taking the Jost function / (k) entering into the well-known S- 
matrix expression [8] 

in the form of a rational function having some simple poles and zeros and exhibiting a correct 
asymptotic behavior at high energies — that is, approaching unity at infinity: 

lim / (k) = 1 . (2) 

fe— >oo 

In the simplest case of only one state in the system, the Jost function has the form 

where the parameter A is always positive, whereas the parameter a is positive for a bound 
state and negative for a virtual state. Upon the substitution of (3) into (1), the S matrix for 
the system having one state assumes the form 

k — la k — i\ 



This state has the energy 



Eo = -^T' (5) 



where m is the reduced mass of the system and U is the Planck constant. The first pole factor 
in expression (4) for the S matrix corresponds to a physical bound or a virtual state of the 
system. The second factor in the S matrix (4) includes the well-known redundant pole [8, 
10], which is associated with no bound state of the system. The redundant pole ensures the 
correct asymptotic behavior (2) of the Jost function. The nuclear vertex constant G 2 for the 
state being considered is directly expressed in terms of the residue of the S matrix at the pole 
k = ia as 

G 2 = m\ 2 Res S (k) , (6) 

k=ia 



where X = Ti/mc is the reduced Compton wavelength of the system. Calculating the residue 
on the basis of (6) and (4), we express the vertex constant in terms of the parameters a and 
A as 

G 2 = -2nX 2 a^±4 ■ (7) 
a — X 

Using (4) and considering that the S matrix is expressed in terms of the phase shift 5 (k) 



as 



S (k) = e 2 ^ = , (8) 

cot — I 

we represent the function k cot S in the form 

kcotS = -- + \r e k 2 . (9) 

This is the expression of the effective-range approximation, with the scattering length a and 
the effective range r e being given by 

11 

a=- + -, (10) 
a A 

2 , N 

a + A 

Thus, the effective-range approximation (9) immediately follows from the Bargmann repre- 
sentation (4) of the S matrix for a system having one state. Formulas (10) and (11) relate the 
parameters of the effective-range approximation to the S'-matrix parameters and, together 
with (5) and (7), the low-energy scattering parameters to the binding energy E and the 
nuclear vertex constant G 2 (parameters of the bound state of the system). 

As a specific example, we will now consider neutron-proton scattering in the triplet spin 
state. In this case, the system features one bound state, the deuteron, with the binding energy 
being = 2.225 MeV. The nuclear vertex constant G 2 d for the deuteron corresponds to the 
d — > n + p vertex and takes the value of G 2 d = 0.43 fm [9]. It should be recalled that the 
nuclear vertex constants are directly related to the asymptotic normalization factors for the 
bound-state wave functions. The latter factors are often introduced in the analysis along with 
the nuclear vertex constants. In the case under consideration, the constant G 2 d is expressed in 
terms of the dimensionless asymptotic normalization factor C d for the deuteron wave function 
as 

G 2 d = 2nX 2 aC 2 d , (12) 



where A" = 2Atv, withA^ = h/m^c being the nucleoli Compton wavelength (mjv is the nucleon 
mass). The relevant numerical value is Cj = 1.673. With the aid of (7) and (12), we can 
easily express the ^-matrix parameter A in terms of the deuteron normalization factor as 

A = 4|±± (13) 



Substituting (13) into (10) and (11), we obtain 



9 r 2 

r e = - ( 1 - ±) . (15) 



a v a 



Formulas (14) and (15) yield explicit expressions for the low-energy np scattering parameters 
in terms of the bound-state parameters, the deuteron binding energy Ed and the deuteron 
nuclear vertex constant G\. For the scattering length and the effective range, the substitution 
of the above experimental values of Ed and G\ into expressions (14) and (15) yields 

a = 5.41/m, (16) 

r e = 1.74/m. (17) 
These values are very close to the experimental values [11] 

a expt = 5.42 fm , (18) 

r e e xpt = 1.76 fm . (19) 

This corresponds to the well-known fact that low-energy neutron-proton scattering can be very 
well interpreted by using the effective-range approximation (9). Thus, low-energy neutron- 
proton scattering in the triplet spin state can be accurately described on the basis of data on 
the deuteron bound state. 
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3. DERIVATION OF THE VAN OERS-SEAGRAVE FORMULA FROM THE 
BARGMANN REPRESENTATION OF THE S MATRIX INVOLVING TWO 
STATES: NEUTRON-DEUTERON SYSTEM 



Let us now consider elastic neutron-deuteron scattering in the doublet spin state at energies 
below the threshold for deuteron breakup. In this case, there are two states in the system, 
the triton ground state and its virtual state. The existence of the latter has been firmly 
established [12]. For the case where the system has two states, we take the Jost function in 
the form of the rational function 

k + iX k + i[i 

which has the correct asymptotic behavior given by (2). The parameters a, A, and /i, which 
appear on the right-hand side of (20), are positive, while the parameter f3, which corresponds 
to the virtual state, is negative. Substituting (20) into (1), we obtain the S matrix for the 
system having two states in the form 

S (k) — ^ + * a ^ + ^ k + i\ k + ifi 
k — ia k — i(3 k — i\ k — ifi 

The energies of the bound and the virtual state are given by 

£ » = -^r- < 22 » 

ft = -^. (23) 

where the reduced mass m is (2/3)m7v- The first and the second pole factor in expression (21) 
for the S matrix correspond to, respectively, the bound and the virtual triton state. The third 
and the fourth factor involve redundant poles of the S matrix. The nuclear vertex constants 
G\ and G\ for, respectively, the bound and the virtual state are directly expressed in terms 
of the residues of the S matrix at the respective poles k = ia and k = if3 as 

G% = m\ 2 Res S (k) , (24) 

k=ia 

G 2 = iixX 2 Res S (k) , (25) 

k=i(3 



where X = §Aat- Calculating the residues and using (24), (25), and (21), we find that the 
vertex constants are expressed in terms of the S-matrix parameters as 



G\ = -2^a ^±4 ^±4 ?±J± , (26) 
a — pa — X a — ji 

<?, = 2^j£±|§±*§±ii. (27) 
a — p p — X p — n 

For a further analysis, it is worthwhile to introduce the following combinations of the 
S'-matrix parameters: 

p = a + p, (28) 

q = a(3, (29) 

u = X + pL, (30) 

v = \fx. (31) 

Multiplying the corresponding factors in the numerator and the denominator of (21), we can 
represent the S matrix in the form 

, . _ k 4 — (q + v + pu) k 2 + qv + ik [(p + u) k 2 — pv — qu] 
k 4 — (q + v + pu) k 2 + qv — ik [(p + u)k 2 — pv — qu] 

Comparing (32) with the representation in (8), we recast the expression for kcotS into the 
form 

qv - (q + v + pu) k 2 + k 4 

kcotd = -— — . (33) 

— pv — qu + [p + u) k z 

Upon dividing the polynomial in the numerator by the polynomial in the denominator, we 
arrive at 

kcot5= -A + Bk 2 - — %— . (34) 

1 + Dk z 

This expression is nothing but the well-known empirical van Oers-Seagrave formula [4], which 
describes well low-energy neutron-deuteron scattering in the doublet spin state. The parame- 
ters in expansion (34) were obtained in [4] by fitting low-energy experimental data. Comparing 
(34) and (33), we express the van Oers-Seagrave parameters in terms of the S'-matrix param- 
eters (28)-(31) as 

A = pq + uv + pu (p + u) ^ 
(p + uf 

£? = — ^, (36) 
p + u 

7 



= qv _ pq + uv+ pu (p + u) 
pt> + (p + M ) 2 

D = --P±^. (38) 
pv + 

Thus, the van Oers-Seagrave formula (34) immediately follows from the Bargmann repre- 
sentation (21) of the S matrix for the system having two states. The same two states in the 
system are responsible for the pole in the expression for the function k cot 5. We note that, in 
the S matrix, it is necessary to take into account, along with two physical poles corresponding 
to the bound and the virtual state, two redundant poles, which ensure the correct asymptotic 
behavior of the Jost function. For the case of potentials leading to a finite number of poles 
in the S matrix, it follows from the Levinson theorem that the number of redundant poles is 
determined by the total number of bound, virtual, and quasistationary states. If we use an 
interaction leading to two states in the system, a ground and a virtual one, the S matrix will 
therefore have two redundant poles. 

Since the neutron-proton system considered above has a single state, a bound state in the 
triplet channel or a virtual one in the singlet channel, the S matrix for a system having one 
state and the corresponding effective-range approximation for the function k cot 5 provide a 
good approximation for the np interaction. The limiting transition from the van Oers-Seagrave 
formula to the effective-range approximation can easily be obtained by recasting formula (33) 
into the form 

, . -l/a + c 2 k 2 + c 4 k 4 

kcot5 = rfm - 2 . (39) 

The two-state S matrix (21) reduces to the one-state S matrix (4) if the second state goes to 
infinity; that is, (5 — > oo and /i — > oo. It can easily be seen that, in this case, the coefficients 
c 4 and D in (39) vanish, so that expression (39) reduces to the effective-range approximation 
(9). Thus, we can see that the van Oers-Seagrave formula is a direct generalization of the 
effective-range approximation to the case of a system having two states. 

The nuclear vertex constants G\ and G\ for, respectively, the bound and the virtual triton 
state are expressed in terms of the corresponding dimensionless asymptotic constants C\ and 
Cy as 

G% = ZidfaCl , (40) 
G 2 V = ZnX 2 {3Cl . (41) 



Using (26), (27), (40), and (41) and taking into account (30) and (31), we obtain 

[ +l) A(a + f3) 2 -Q(a + (3) 2 (C 2 +C 2 ) + 9(a-(3) 2 C 2 C 2 ' 1 ' 

4 (q + (3) 2 + 6 (q + l) 2 (Cg. + Cg) + 9 (a - /?) 2 Cfg 
aP 4 (q + /3) 2 - 6 (q + (3) 2 (C| + C„ 2 ) + 9 (q - /3) 2 C%C* ' 
Formulas (35)-(38), together with (28), (29), (42), and (43), provide explicit expressions for 

the van Oers-Seagrave parameters in terms of the parameters of the bound and the virtual 

triton state (their energies and nuclear vertex constants). Since these general expressions are 

very cumbersome, we consider here only the expressions for the scattering length and for the 

pole of the function k cot 5. 

4. EXPRESSIONS FOR THE nd SCATTERING LENGTH AND FOR THE 
POLE OF THE FUNCTION kcottf IN TERMS OF THE PARAMETERS OF 
THE BOUND AND THE VIRTUAL TRITON STATE 

As can be seen from (34), the doublet nd scattering length a = 2 a n( i is expressed in terms 
of the van Oers-Seagrave parameters as 

° = ~A~hC' < 44 > 

or in terms of the S'-matrix parameters as 

P u 1 111 
a = Z + - = - + - + - + -. 45 

q v a p A fi 

Substituting (28), (29), (42), and (43) into (45), we obtain 

a = 6 (l + i\ 2 ( Q ± g) (ggj ± aC %) ~l~ 3 (q P) 2 C^Cl 

\a p) 2(a + (3) 2 [2 + 3(C*+C2)]+9(a-(3) 2 C*Cf 1 ' 

This formula expresses the doublet nd scattering length in terms of the parameters of the 

bound and the virtual triton state. From (38), we similarly find that the pole k% = l/D of 

the function k cot 5 is given by 

k 2 = a 2 (q + g) (gg| + aCp + 3 (q - P) 2 C 2 T C 2 V 
P 2(a + (3)(aC 2 +(3C 2 )-3(a-/3) 2 C 2 Cf 1 1 

We note that the effective range for nd scattering can be expressed in terms of the van 

Oers-Seagrave parameters as 

r e = 2(B + CD) . (48) 
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Along with the effective range, which is anomalously large, r e ~ 500 fm, for doublet nd 
scattering, the so-called amplitude slope parameter [6], given by the dimensionless quantity 
|a^a 2 r e with wave number a d = 0.2316 fmT 1 corresponding to the deuteron, is often used in 
the literature. 

5. NUMERICAL CALCULATION OF THE LOW-ENERGY nd 
SCATTERING PARAMETERS AND THE CONSTANT C 2 V 

With the parameters specifying the bound and the virtual triton state, we have calculated 
numerically the low-energy nd scattering parameters a, Eq = — and -^a d a r e by formulas 
(46)-(48). The calculation has revealed that these parameters of nd scattering are weakly 
sensitive to variations in the energies E and E v of the ground and the virtual state, respec- 
tively, and to the asymptotic normalization factor C\ for the bound state. At the same time, 
these parameters greatly depend on the asymptotic normalization factor C% for the virtual 
state. The results of the calculations for the low-energy nd scattering parameters are quoted 
in the table for various values of the constant C%. In the calculations, the parameters of the 
triton were set to the values of E T = 8.48 MeV [13] (E T = \E \ +e d ), B v = \E V \ = 0.482 MeV 
[12], and = 3.5 [1, 2, 6], which follow from a direct analysis of experimental data. 

From the results quoted in the table, it can be seen that the parameters of nd scattering are 
highly sensitive to variations in the constant C 2 . The experimental value of C% = 0.0504 [12] 
adopted at present leads to the nd scattering length a = 1.03 fm, which differs considerably 
from its experimental value now established to a high precision. The experimental value of 
the nd scattering length [14], 

a expt = 0.65 fm, (49) 
can be fitted with the asymptotic-constant value 

Cl = 0.0592 , (50) 

which agrees well with the value of Cl = 0.06 calculated in the two-body model of nd in- 
teraction simulated by the Hulthen potential [2]. Thus, the numerical value in (50) for the 
asymptotic normalization factor C% for the virtual state must be treated as the result of the 
theoretical calculation on the basis of our model. The high sensitivity of the scattering param- 
eters to this quantity and a deviation of the theoretical results from its available experimental 
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value indicate that a refinement of the experimental value of C 2 is necessary. The above 
highlights the crucial importance of the factor C 2 for exploring the properties of nd scattering 
in the doublet spin state. The experimental value in (49) for the nd scattering length and, 
accordingly, the value in (50) for C 2 lead to the following values for the pole of the function 
k cot 5 and the amplitude slope parameter: 

£ = -0.1313 MeV, (51) 

l -a\a 2 r e = 1.52 . (52) 

The value calculated for the pole position e and presented in (51) agrees well with e = 
—0.15 MeV quoted in [15] as an experimental value. The slope parameter (52) complies fairly 
well with the value of \a\a 2 r e = 1.35 [6] calculated theoretically from the Faddeev equations 
with separable nucleon-nucleon potentials. 

For the van Oers-Seagrave parameters, a numerical calculation with C 2 from (50) yields 



A = 0.3198 /m" 1 , (53) 

5 = 0.7698/771, (54) 

C = 1.2190 fmT 1 , (55) 

D = 236.927 /m 2 . (56) 



The function k cot 5 calculated in the van Oers-Seagrave approximation with the parameter 
values (53)-(56) is shown in the figure versus the energy k 2 . It can be seen that the theoretical 
curve faithfully reproduces experimental data from [4]. In summary, we have obtained a 
complete description of low-energy doublet nd scattering in terms of the parameters of the 
bound and the virtual triton state and demonstrated that this description is consistent with 
experimental data. 

6. CONCLUSION 

It has been shown that the van Oers-Seagrave formula for doublet nd scattering — it 
was originally deduced from a purely empirical consideration — immediately follows from 
the Bargmann ^-matrix representation corresponding to the presence of two triton states, a 
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bound and a virtual one, in the system. That the pole term proves to be necessary in the 
function k cot S is also due to the presence of two states in the system. Our analysis has given 
simple expressions relating the van Oers-Seagrave parameters to the S'-matrix parameters. It 
has been shown that the van Oers-Seagrave formula is a direct generalization of the effective- 
range approximation to the case where the system being considered has two states; when one 
of these states goes to infinity, the former reduces to the latter. 

For the doublet nd scattering length and the pole of the function k cot 5, we have obtained 
simple explicit expressions in terms of the energies of the bound and the virtual triton state 
and the nuclear vertex constants for these states. The resulting formulas make it possible 
to perform numerical calculations and to analyze the parameters of low-energy nd scattering 
versus the parameters of the triton. In particular, the calculations have revealed that the 
parameters of nd scattering are highly sensitive to variations in the asymptotic normalization 
factor C% for the virtual triton state. The experimental value of the nd scattering length is 
fitted at the value of C 2 U = 0.0592, which differs from the experimental value presently accepted 
for this constant. Therefore, the constant Cl is of paramount importance for studying the 
properties of the nd system and needs further experimental refinement. 
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Function k cot 5 calculated for doublet nd scattering in the van Oers-Seagrave approximation 
with the parameter values (53)-(56) versus the energy k 2 . The experimental data were taken 
from [4]. 
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